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Integration of a System of Ordinary Differential Equations' 

up to a Specified Value of One Variable. 

Closed with 3 program parameters . 

60 

0, 1, 2 

This subroutine will integrate a system of n simultaneous 

first order ordinary differential equations, in which each 

derivative is expressed explicitly in terms of the variable* 

thus: 

y = f ( y o ,y l* ""* y n-l^ 

y, i = f i fr ' y i> •••' y B .i) 



^'n-l f n-l ( y o' y i> 



METHOD OF USE 



"' y n-l> 



Integration is carried out by means of the Gill-Kuita method, 
described in routine Fl from. a set of initial values of the 
variables, until one specified variable y fc reaches the value 
Y^. Equal step intervals are used until y_ is near to Y . 
At this tijae the step interval h is automatically readjusted 
in accordance with the rule *=(\-y±)/y\ of the Kewton- 
Raphson method until y k converges to Y fc with an error of 
at most 2" . Control is then returned to the main routine. 
Locations 3 through 7 must contain the preset parameters a, 
b, c, n, 4 as integers before and during input of routine F$. 
We must have a < b < c . 

a is the location for the first variable y . other 
variables occupy succeeding locations so that y 
occupies a + i. The initial values are originally 
placed in these locations and routine F5 alters them 
during the integration process . 
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b is the lnatior. for 4 :h.e .first sealed «?.'.' vat Its 

2 m h y f - Other derivatives occupy succ~?diw^ 
locations sc that 2 h y. ' occupies b + i yy- 
auxiliary subroutine must calculate these i ' 
v&tives and place them in their locations , 

c Locations c + i, (i = 0, 1, ■■ ■■ . s n-l) are used 

as temporary storage by F9» These locations 
must be cleared before entering F$ with a new 
set of initial values . 

n is the number of differential equations to be 
solved . 

d is the location of the first word of the auxiliary 
subroutine . 

Entry to F*> is made by means of the instructions 



p 


50 mF 




50 p 


p+1 


26 - 




00 (a+k)F 



and control is returned to the left hand side of p + 2. Tfc.e 

scaling parameter m Is used in the same fashion as in routine 

Fl and, if varied, will change the step interval h provided 

the auxiliary subroutine is unchanged . The program parameter 

a + k is the location of y, , the variable which is brought 

to the value Y,_. The value Y, must be in the accumulator at 
k ■£ 

the time of entry. 

The auxiliary subroutine must be provided by the user. It 
must be programmed as a closed subroutine using standard 
entry and must use the variables obtained from a,a+l, • .»,a+n-l 
to compute each 2 m h f and place it in location b + i, for 
each i = 0,1, . ..,n-l. 
HQTE8 A more complete description of the method of integration is 

contained in the description of routine Fl. For a discussion 
of accuracy, scaling, and the choice of parameters h and m 
the reader is referred to this description. 



EXAMPLE 
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The time per step is approximately the same as for Fl but 
the number of steps per entry into 5$ may be difficult t.-> 
determine. In case y fc never reaches Y, the integration 
process will not terminate. 
One wishes to find values of x such that 

.2 



, 'x -t 



VST 



dt 



assumes the values .9, .99, and .999. First write the 
integral as a set of first order differential equations by 
introducing a new variable 



u = 



2 



We may then write 

x' = 1 

u* = -xu 

y 1 = u s 

as the set of equations to be solved by 1$ ,_"/' Initial conditions 

x = 0, u = — = , y = - may be placed in S3, 1S3, 2S3 before 
y2jc 

starting the integration. The entry 



p 


50 mF 
50 p 


p+1 


26 - 
00 S3 



may be made repeatedly with -.9, .99, and .999 in the accumulator 
to obtain the corresponding three values of x. The initial 
conditions need not be reset.. 
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26 1000N (F5) 

h 







L5 5F 
IA 6f 






1 


00 20F 


1 


-cfn in 53 




^6 53L ' ! 


i 


2' 


L5 3F j fj 




k2 5kL | j 


• 


3 


00 20F j j 


5-a in 5^ 


k 


> k6 5k L | j_ 
L5 ^F ; j 

! LO 3F ; ! 

f ' i 


1 
i 


5 
j 


! k2 55L ' j 
| 00 20F ' j 


-b-a in 55 

1 
i 


6 


j k6 55L ; I 






1 L5 5F i # r 


■^ 


7 


LO kF ; 






U2 56L 1 


-c-b in 56 


8 


00 20F ; I 




9 


k6 56L ' ! 
1 L5 58L ; 






! 22 IOIVf! 


End of interlude 


10 


kO IF 




1 


; LO 59L ! 


y k m 1 


! ii 


; kO 2F 




i 


L3 F j 


V Y k in 2 


i 12 


j 10 (m)F 




|. 


| L6 2F 


is |Y k -y k | < h y k 


! 13 

i 

i 


! 32 38L 
L3 2F 


is |Y k -y k l < 2" 59 ? 


ik 


F^ 5TL ' 

: 36 (p+2)f 




15 


( 

50 57L 
L5 2F 




r 

1 

i 

1 


; 16 


00 (m)F 
66 F 




I 

s-K - y t 

j! hy k 




' 


. ■.. . , 


*.^_ — • ■ - - - - 
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.. 17 


S5 F 








•■ 1 ■ ■ ■ 


^0 58L 








18 


22 38L 
50 58L 








*. 19 


L5 (q, JF 

^0 IF 








20 


79 ( kijJ )F 

40 F 








21 


LO (1)F 








22 


LO IF 
kO 2F 
50 (A,)F 

o 




(k l,J" B J 4 W ) 




23 


U 2F 
L4 2F 








2k 


10 (m)F 
kO 2F 




f k U " b j "kjJCAj+D.a- 


i 


25 


>M yij )F | 






26 


40 (y i,J + l 

L5 2F 
50 57L 


)F 


step y ± 




27 


00 (m)F 
kO 2F 




fTOmT l s t + l 




28 


50 (c )F 

7J F 








29 


LO F 








30 


Lk IF 
Lk 2F 
Lk 2F 




C . k, . ■ - k . . 




31 


L*J- 2F 










40 ( *i, J+ i 


IF 


q iJ + 3 r i,j+l 




32 


L5 25L 
L4 29L 








33 


k2 25L 
^6 25L 








3^ 


L4 55L 
k6 20L 









F5 



LOCATION 
35 

56 

37 
38 
39 
40 
41 
42 
43 

45 
46 

47' 

48 

49 

50 

51 



ORDER 

L4 56L 

42 31L 

46 19L 

LO 53L 

32 18L 

50 37L 

26 S7 

L5 (49 )L 

42 38L 

32 8L 

46 21L 

10 10F 

L4 34l 

42 22L 

46 28L 

L5 54L 

26 33L 

00 F 

40 F 

00 F 1/2 

NO F 

00 F -1/2 

40 F 

oo 207106781186J 

80 F 

00 292893218814J -1/ fi~~ 

80 F 

00 1666666S6667J -5/6 ' 
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•Step addresses depending on i until i = n 



F5 



-Increase j and adjust addresses depending 
on j 



A, 







LJ 1025F 
06"1074l 
LJ 3129L 
06 3123L 
LF 57L 
06 2100L 



-11, 1, 

25, 50L 

-9,2,27, 
51L 

-8,2,26, 
52L 



te 



C 1 , Ag 



A l> C 2 



„._JL 






LOCATION 

53 

5 4 

55 
56 
57 
58 



3 

5 
6 

7 
8 

9 



L 



ORDER 

LJ 1025F 
07 53L 
64 (c+n)F 
00 49L 
80 (a)F 
00 (a)F 
00 (b-a)F 
00 (b-a)F 
00 (c-b)F 
00 (c-b)F 
00 F 
00 F 
00 F 
26 L 

26 IN 

40 59L 
K 5F 
46 12L 
46 i6l 
46 24l 
46 27L 
42 5L 
F4 57L 
4a. 14L 
LJ 44l 
40 58l 
L5 (p+l)F 
42 9L 
L4 55L 
42 8L 
50 7L 
26 S7 
L5 (b+k)F 
40 F 

L5 (a+k)F 
01 60K 
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F5 



-11,1,28, 
53L 



zero 



End constants set by interlude 



Y k in 59 



-m in 12, l6, 24 and 27 

p+1 in 5 
-link in 14 

-K in 5.8 

-a+k'in 9 

~-b+k in 8 

-Enter Auxiliary 



|-2 m h y* k in 



